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Summary. We report on large-scale Monte Carlo simulations of the three-dimen- 
sional 4-state Potts ferromagnet subject to quenched, random bond dilution. For 
small dilutions the rather strong first-order phase transition of the pure system is 
found to persist, whereas for larger dilutions the theoretically expected softening 
to a continuous transition is confirmed. The properties of the underlying disorder 
distributions of thermal observables are discussed and illustrated with a few selected 
examples. 



1.1 Introduction 

Physical systems with quenched, random impurities display under certain 
conditions a completely different behaviour than the pure systems. One (of- 
ten unwanted) experimental realization are point defects which may be mod- 
elled theoretically by random site dilution [1]. If the pure system exhibits a 
continuous phase transition, the influence of quenched, random disorder may 
drive the system into a new universality class, provided the critical expo- 
nent a p of the pure system's specific heat is positive (Harris criterion) [2]. 
Also in the case of a first-order phase transition in the pure system, quenched 
disorder can have a dramatic effect: Phcnomcnological renormalization group 
arguments suggest the possibility of a softening to a continuous transition [3] . 
While in two dimensions this effect has been theoretically proven [4] and nu- 
merically confirmed [5] for any small amount of disorder, in three dimensions 
(3D) one expects that with increasing disorder the strength of the transition 
is gradually weakened until beyond a (tricritical) concentration of impurities 
a real softening to a second-order transition sets in [6,7]. 

Recent Monte Carlo (MC) computer simulations of the 3D siie-diluted 
3-state [8] and frond-diluted 4-state [9] Potts models, which both exhibit 
first-order phase transitions in the pure case, have confirmed this expecta- 
tion numerically. For a complementary study using high-temperature series 
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Fig. 1.1. Distribution of the magnetic susceptibility X{J} f° r P — 0-56 and L = 96, 
as obtained at a temperature close to the point where the average [X{j}(T)]a,v has 
its maximum. The solid line in (a) is the running average, and the vertical lines in 
(b) show from left to right the most probable, median, and mean value. 



expansions, see Ref. [10]. One important issue in MC studies is the number 
of independent disorder realizations that are required for reliable quenched 
averages. An answer to this question can usually only be given a posteriori, 
after having already estimated the associated distributions of thermodynamic 
quantities such as the susceptibility, since these distributions can be highly 
asymmetric. As is illustrated in Fig. 1.1, one result of our large-scale simula- 
tions [9] is that this asymmetry is particularly pronounced in the softening 
regime. In what follows we shall further illustrate and elucidate this effect 
and its consequences with a few selected examples. 



1.2 Model and Simulation Setup 

The 3D bond-diluted g-state Potts model is defined by the Hamiltonian 

H = -^JijSai^/, 0"i = l, (i-i) 

(V) 

where the sum extends over all pairs of neighbouring sites on a cubic lattice of 
size L 3 with periodic boundary conditions, and the couplings Jij are chosen 
according to the distribution 

piJij) = p6(J i:i - J) + (1 - p) S^), (1.2) 

where p is the concentration of magnetic bonds such that p — 1 corresponds 
to the pure case where the 4-state model exhibits a rather strong first-order 
transition at K c (l) = K c = J/{k B T c ) = 0.62863(2). 

We simulated the diluted model in the regime p c < p < pt of second-order 
transitions, where p c — 0.248 812 6(5) is the bond-percolation threshold [11] 
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Fig. 1.2. Temperature variation of four different realizations at p — 0.56 and 
L = 40 (long and short dashed lines). The solid line with a maximum around 
X{J} ~ 160 shows the average susceptibility \ = [X{./}]av- 



and p t « 0.80 denotes the tricritical concentration [9], with the Swendscn- 
Wang cluster algorithm [12], and in the regime pt < P < 1 of first-order 
transitions with the multibondic method [13]. The numerically determined 
phase diagram in the dilution-temperature plane was found to agree very well 
with the single-bond effective-medium (EM) approximation [14], 



(l-Pc)e^-(l-p) 



P-Pc 



(1.3) 



where K c (l) is the transition point of the pure system and p c the percolation 
threshold given above. Our analysis of both the autocorrelation time and 
the interface tension [9] led to the conclusion of a tricritical point around 
p = {)M. 

To arrive at these results, for each dilution, temperature and lattice size, 
the MC estimates {Q{j}) of thermodynamic quantities Q{jy for a given ran- 
dom distribution {J} of diluted bonds were averaged over 2 000 - 5 000 dis- 
order realizations, 

Q = [(Q{j})U = Jvj t3 p{J l3 )(Q { . J} ) = jd(Q {J} )V((Q { j } ))(Q { j } ), (1.4) 

where V((Q{j})) denotes the empirically determined distribution of {Q{j}) 
discussed in the next section. 



1.3 Distribution over Disorder Realizations 

In Fig. 1.1 the distribution of susceptibility measurements at a temperature 
close to the peak of its disorder average is shown for p = 0.56 and L = 96, i.e., 
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Fig. 1.3. Time evolution of the magnetization for measurements close to the 
(average) susceptibility peak for two different realizations at p = 0.56 and L = 40, 
and the associated thermal equilibrium distributions. 

for a dilution in the centre of the softening regime 0.25 w p c < p < p t ~ 0.80. 
The solid line Fig. 1.1(a) is the running average which, despite the quite 
asymmetric distribution, appears rather stable from about 400 realizations 
on. The resulting probability density (histogram) and distribution (accumu- 
lated density) are depicted in Fig. 1.1(b). By the vertical solid and dashed 
line we have also indicated the average and median value, respectively. The 
dotted vertical line at about X{J) = 350, which perfectly coincides with the 
most probable value, has been computed by averaging over only those 50% 
of the realizations whose X{j}- v alue is smaller than the median value. 

In Fig. 1.2 wc illustrate how the average susceptibility emerges from the 
individual contributions of each disorder realization. The location of its max- 
imum is indicated by the vertical line. Two events corresponding to large 
values of X{J} (rare events) are shown as long dashed lines, and two events 
with small values of X{j} (typical events) are shown as short dashed lines. 

Since each X{j} reflects the thermal fluctuations in the magnetization m, 
we also had a closer look at the time evolution of m and the associated equi- 
librium distribution. Two examples are shown in Fig. 1.3, where the upper 
part corresponds to a rare event with large susceptibility and the lower one to 
a typical event with smaller magnetization fluctuations. The first realization 
still behaves first-order like as can be seen more clearly on the right part of the 
figure, where the resulting thermal probability density of the magnetization 
exhibits a double peak. 
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1.4 Conclusions 



By performing large-scale Monte Carlo simulations we have investigated the 
influence of bond dilution on the thermodynamic properties of the 3D 4-state 
Potts model. For strong dilution with p < 0.80 we obtained clear evidence 
for softening to a continuous transition. Here we have focused on a rather 
detailed study of the thermalization properties of individual realizations and 
the disorder distribution of thermodynamic quantities such as the suscep- 
tibility. Among other observations this reveals that in the softening regime 
a few rare realizations may still exhibit signals reminiscent of a first-order 
transition, e.g., a double peak in the magnetization histogram. On the one 
hand, this gives rise to rare events with large susceptibilities leading to quite 
asymmetric distributions over the disorder, and on the other hand this com- 
plicates the technical issue of deciding which Monte Carlo update algorithm 
is best suited for the problem at hand. 
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